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Aims. Large HI shells, with diameters of hundreds of pc and expansion velocities of 10-20 kms 1 have been detected in their hundreds in 
f**-) , the Milky Way and are well observed features of local gas rich galaxies. These shells could well be predicted as a result of the impact of OB 
associations on the ISM, but doubt has been cast on this scenario by the apparent absence of OB stars close to the centres of a large fraction 
of these shells in recent observations of the SMC. Here we present observational evidence within a energetically consistent framework which 
strongly supports the scenario in which OB associations do produce the giant HI shells. 
O ' Methods. Using Fabry-Perot scanned Ha emission line mapping of nearby galaxy discs, we have detected, in all the H n regions where the 
observations yield sufficient angular resolution and S:N ratio, dominant Ha shells with radii a few tens of pc, expanding at velocities of 
t/\ • 50-100 kms -1 , and with gas masses of 10 4 - 10 5 M o . In previous studies, we found that stellar winds alone can account for the energetics of 
most of the Ha shells, which form initially before the stars explode as SNe. We have applied a simple dynamically consistent framework in 
which we can extrapolate the properties of the observed Ha shells to a few 10 7 yr after the formation of the OB stars. The framework includes 
the dynamical inputs of both winds and SNe on the surrounding ISM. The results give quantitative statistical support to the hypothesis that the 
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Ha emitting shells are generic progenitors of the HI shells. 
^ Results. The results are in good agreement with the ranges of masses (~10 6 M o ), velocities (up to ~20 kms -1 ), and diameters (up to ~500 
' pc) of representative HI shells observed in nearby galaxies. The combined effects of stellar winds, acting during the first few 10 6 yr, and SN 
explosions, "switching on" subsequently, are required to yield the observed parameters. 
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1. Introduction 

The presence of neutral hydrogen shells in the gas of disc 
galaxies is a widely known and observed phenomenon. Starting 
over 25 years ago Heiles (1979, 1984, 1990) detected large 
shell structures in Galactic neutral hydrogen, and in the inter- 
vening period a number of authors (e.g. Brinks 1981, Brinks 
& Bajaja 1986, Deul & den Hartog 1990, Puche et al. 1992, 
Kamphuis & Sancisi 1993, Stavely-Smith et al. 1997, Wilcox 
& Miller 1998, Walter & Brinks 1999, Kim et al. 1999, Wills 
et al. 2002) have shown that such shells, often detected around 
holes in HI, are present in nearby gas rich galaxies. More re- 
cently Mclure-Griffiths et al. (2002) mapped hundreds of HI 
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shells within the Galaxy, and Ehlerova & Palous (2005) showed 
that in the outer galaxy these shells occupy a minimum of 5% 
of the disc volume. Although telescope resolution limits the 
minimum detectable size of an HI shell to some tens of pc even 
in nearby galaxies, many have been detected on scales of hun- 
dreds of pc, up to 1 kpc in diameter. 

Brinks & Bajaja showed that in M31 HI holes smaller than 
300 pc are correlated with OB associations, but no reliable cor- 
relations could be found for bigger holes. However, when suf- 
ficient spatial resolution can be applied, such as in the case of 
the SMC, it is notable that positional coincidence of the HI 
shell with a corresponding OB association occurs in only a mi- 
nority of cases, as shown by Hatzidimitriou et al. (2005). These 
authors showed that of the 509 expanding HI shells catalogued 
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in the SMC only some 40% are spatially correlated with OB 
associations, and some 10% lie in low density fields with no 
associated young stellar objects. Rhode et al. (1999) found evi- 
dence of young stellar clusters in only a small fraction of the HI 
holes in the dwarf galaxy Ho II. The age of the central stellar 
cluster located within the supergiant shell in the dwarf irregu- 
lar galaxy IC 2574 is in agreement with the kinematic age of 
the shell, supporting the scenario in which the supernovae and 
stellar winds create the giant HI holes (Stewart & Walter 2000). 

Mechanisms of different kinds for producing large HI shells 
have been presented in the literature. The possibility that the in- 
fall of a high velocity cloud can give rise to HI holes and shells 
has been discussed by a number of authors with some quite de- 
tailed models (e.g. Tenorio-Tagle & Bodenheimer 1988, Rand 
& Stone 1996, Santillan et al. 1999, Murray & Lin 2004). In the 
last two of these studies magnetic fields were also included. 
Another possibility, proposed by Hatzidimitriou et al. (2005) 
is a turbulent origin for at least some of these structures. The 
present article does not claim the scope to contrast those the- 
oretical models with observations. All we will do is to present 
our scenario, based on observations of expanding Ho- shells in 
H ii regions, which offer a new quantitative element capable of 
contributing to explain the range of parameters observed for HI 
shells, including the largest of them. 

The essential point of this article is to present evidence, 
using a set of observational data not hitherto applicable to 
this scenario (that of the presence and properties of expand- 
ing ionized shells within luminous H n regions) within an en- 
ergetic framework of a set of models, that these ionized shells 
are the precursors of the previously observed class of neutral 
shells which have larger radii and lower expansion velocities. 
Although the scenario we present here is conventional in the 
sense that it links the HI shells to OB associations, our quan- 
titative treatment is useful in showing that observations such 
as those of Hatzidimitriou et al. (2005) are in fact compatible 
with it. Typical measured velocities of HI shell expansion (e.g. 
Walter & Brinks 1999), are of order 10 kms -1 , much lower 
than the velocities detected in OB stellar winds (Herrero et al. 
1992) and very much lower than velocities in supernova rem- 
nants (Weiler & Sramek 1988, Chevalier 1977). The expansion 
velocities we find for the Ho- shells are of order of 100 kms 1 . 
In the context of our schematic model developed in this arti- 
cle these velocities are characteristic of the successive stages 
of the developing shell as it expands and does work against its 
surroundings. The timescales implied for a shell to reach radii 
of order 1 kpc are at minimum a few times 10 7 yr, which gives 
time for the stars in a coeval OB association to disappear. We 
also note that an expansion occurring in an anisotropic zone of 
the interstellar medium (ISM) will expand preferentially in the 
low density direction, so that an initiating star cluster would not 
lie close to its centre; for the largest shells this effect would be 
enhanced by Galactic differential rotation. 

For these reasons we set out to detect the counterparts of 
HI shells at a much earlier stage, while they are more com- 
pact but expanding at higher velocities. Of course this implies 
that we cannot be looking at the same objects as those detected 
via the HI shells and holes, but rather at the H n regions which 
surround the OB associations and which mark the energetic im- 



pact of these associations on the ISM. To do this we used the 
optical technique most akin to 21 cm line mapping of HI: two 
dimensional emission line mapping in Ha with a Fabry-Perot 
interferometer. Some of these results have been reported pre- 
viously (Relano & Beckman, 2005a, Paper I), and showed the 
presence, within luminous H n regions, of dominant expansive 
components whose properties indicated that they could be the 
precursors of large HI shells. In section 2 of the present article 
we summarize those observations, and describe further work of 
the same basic type on H n regions in the nearby galaxy NGC 
5194 (M51). In section 3 we describe a very simple dynamical 
model in which the inputs from stellar winds and from super- 
novae are coupled to the surrounding gas, taking energy and 
momentum conservation as basic constraints. We apply this to 
our observed HII regions, using the measured Ha luminosi- 
ties to estimate the massive stellar content used to calculate the 
energy and momentum inputs to the shells. We go on to com- 
pare these results with the range of energies and momenta in 
HI shells from the literature. In section 4 we discuss general 
physical considerations relevant to modelling the implied evo- 
lutionary sequence and give our conclusions. 

2. Evidence for expanding shells in luminous HII 
regions 

Using galaxies whose populations of Hn regions have pub- 
lished photometric catalogues in Ha (Rozas et al. 1996; Rozas 
et al. 2000; Relano et al. 2005b) we used complete intensity- 
velocity maps, also in Ha, made with the TAURUS-II Fabry- 
Perot system on the 4.2m WHT (Paper I) to make an exhaustive 
study of the integrated emission line profiles of the H n region 
population in three spiral galaxies: NGC 6951, NGC 3339 and 
NGC 1530. We found that more than one third of the regions 
in each galaxy showed high velocity low intensity peaks in the 
wings of the main emission lines. The fraction of the line pro- 
files showing these wing features increased with the S:N of the 
observations, so that, for example, NGC 1530, the galaxy with 
the highest S:N was the galaxy with the highest fraction of its 
profiles exhibiting wing features. 

The interpretation we gave (Paper I) to these low intensity 
high velocity features is that they are evidence for supersonic 
expanding shells within the H n regions, and we could use their 
observational parameters to derive the properties of the shells. 
Fitting the features with Gaussians we could derive, first, the 
emission measure corresponding to the feature, which was used 
to derive the column density of a given shell and, using reason- 
able assumptions, go on to estimate its density (see Eq. (1) of 
Paper I). Secondly, the mean velocity separation between the 
central Gaussian emission component and the lower intensity 
wing features gave us the velocity of the expanding shell. 

Here we have taken the dynamical parameters (mass and 
kinetic energy) of a sample of H n regions from Paper I. The 
sample was selected to cover the full range of Ha luminosities 
of the Hn regions observed in the three galaxies (see TableHJ. 
The estimated shell radii shown in TableQ]are rather crude esti- 
mates, since at the distances of the galaxies the linear resolution 
does not allow a clear measurement. Based on previous studies 
of expanding shells in extragalactic luminous H n regions much 
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closer to the observer, 30 Doradus (Chu & Kennicutt 1994) 
and NGC 604 (Yang et al. 1996), in Paper I we used a canon- 
ical shell radius of R s h e ii = 0.2 R reg (i.e. 20% of the Hn region 
radius) and a shell thickness of 4.5 pc. From these estimates 
we could go on to estimate the shell mass as a fraction of the 
total gas mass in the H n region. Average values for the H n re- 
gions in the three galaxies ranged from 20% for NGC 1530 to 
26% for NGC 3359, i.e. of order 20-25%. We could then use 
the estimated mass and the measured velocity to derive an esti- 
mated kinetic energy of expansion for each shell. We obtained 
the equivalent number of 03(V) stars using the Ha luminosity 
of a given region and the value given by Vacca et al. (1996) for 
the ionizing photon output of an 03(V) star. We can later use 
the equivalent number to estimate the kinetic energy input from 
the stellar winds, using the estimate by Leitherer (1998) of the 
wind luminosity for an 03(V) star, integrated over the time the 
star is on the main sequence. This input energy is listed in the 
last column of Table [TJ It is clear that we could have added 
the extra sophistication of looking at an IMF, instead of choos- 
ing a specific stellar spectral type, but the uncertainties implied 
would not have made this option a better one in terms of accu- 
rate quantification of the phenomena. 

2.1. Observations of NGC 5194 

We have followed the method used in Paper I and outlined in 
the paragraph above to study a sample of H n regions in the 
nearby galaxy NGC 5194 (M51). The Ha Fabry-Perot obser- 
vations for this galaxy were obtained on the 1.6 m telescope of 
the Observatoire du Mont Megantic (OMM) using the instru- 
ment FaNTOmM (Hernandez et al. 2003; Gach et al. 2002). 
The final data cube consists of 48 spectral channels with a cor- 
responding total integration time per channel of 5.16 min. The 
channel resolution is 0.15 A (corresponding to 6.94km s~') and 
the effective resolution of the observations is R= 18609, which 
corresponds to 16 kms" 1 . The central wavelength of the red- 
shifted narrow band blocking filter for Ha is 658 1 A, with a half 
width of 9.9A corresponding to 450 km s _1 . 

The Photon Counting Camera had a scale of 1.6"/pixel and 
the seeing during the observations was 3-4". Although this see- 
ing cannot be considered as very good, the high spectral res- 
olution allows us to study well the high velocity features in 
apertures commensurate with the seeing, which allows us to 
study the largest and most luminous H n regions adequately. A 
detailed description of the observations of NGC 5194 is given 
in Daigle et al. (2006a) and the data reduction is explained in 
more detail in Daigle et al. (2006b). The Ha luminosities for 
the Hn regions in NGC 5194 were taken from a catalogue ob- 
tained by Rand (1992), and kindly supplied by the author. The 
radius of a region found in the catalog was used to estimate the 
highest aperture for the line profile of any given region under 
study. 

In Fig. Q] we show two examples of line profiles for Hn 
region 403 and 312 of M51 in Rand's catalogue with the cor- 
responding Gaussian components fitted to this spectrum. These 
figures clearly show the typical red and blue components called 
here wing features. For an extended sample of the line profiles 



from H ii regions in other galaxies the reader is referred to Paper 
I. The results of the analysis for the emission line profiles of the 
H n region sample in NGC 5 194 are shown in TableQ] For each 
H n region we extracted line profiles with apertures of increas- 
ing radius, and fitted them with the corresponding Gaussian 
components. The profile whose wing features yielded the best 
S:N was taken as the representative spectrum of the H n region. 
The observational parameters of the wing features for each H n 
region do not vary significantly when increasing the aperture, 
which allows us to take the best S:N line profile as the rep- 
resentative one for each H n region. The corresponding obser- 
vational parameters, the Emission Measure and the velocities 
of the wing feature Gaussians with respect to the central emis- 
sion peak were used to obtain the dynamical parameters of the 
shell for each region using the same criteria as those in Paper 
I. These parameters are listed in Table Q] The relative errors 
for the observational parameters are 0.2-0.3 for the Emission 
Measures and 0.05-0.08 for the velocity separation between the 
central and the shifted components. It is interesting to note that 
the absolute errors in the fitting procedure for the velocity of 
the shifted components are 2-3 kms 1 , much lower than the 
spectral resolution of the observations. This gives an idea of 
the validity of the Gaussian decomposition procedure. 

3. Energetics of the expanding shells 

3.1. Starburst99 Models 

Our aim here is to extrapolate from the properties of the shells 
we have observed within the Hn regions to later times, tak- 
ing the shells as representative, to show how their parameters 
would evolve with time and then to compare the extrapolated 
shells with HI shells observed in the Galaxy and nearby galax- 
ies. In order to infer the energy inputs from the stars in the 
form of both stellar winds and supernovae, we have used from 
the literature the Starburst99 models (Leitherer et al. 1999). We 
take the case of an instantaneous starburst with a fixed stellar 
mass and integrate using two different Initial Mass Function 
(IMF) exponents, 1.3 and 2.3 for mass intervals (0.1-0.5)M Q 
and (0.5-100)M Q , respectively. We take the conventional as- 
sumption that 8M is the lower limit of a Zero Age Main 
Sequence star which will eventually explode as a supernova, 
and use the Padova evolutionary tracks at solar metallicity. 
Starburst99 then gives, from an initial time of 10 4 yr the (loga- 
rithmic) energy (in erg s _1 ) input due to both stellar winds and 
supernovae as a function of time. 

The principal input parameter in our Starburst99 models is 
the stellar mass associated with the burst, which was inferred 
from the Ha luminosity of the H n region, as described in sec- 
tion 7.2.1 of Relano et al. (2005b). We use the Salpeter IMF to 
estimate the ionizing mass of the star cluster and compare it to 
the mass estimate using an equivalent number of 05(V) stars. 
The normalizing factor obtained is then used when integrat- 
ing the IMF within the appropriate mass limits (0. 1 M - 1 00M o ) 
and the total stellar mass in the cluster is derived. The output 
of Starburst99 gives the energy produced by the stellar winds 
and supernovae within the star cluster as a function of the age 
of the cluster. We have taken the total energy produced by the 
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Table 1. Shell parameters of a subsample of Hn regions in NGC 1530, NGC 3359, NGC 6951 and NGC 5194. The shell 
parameters of the Hn regions in the first three galaxies are from Paper I, while the rest are derived here. Column 1: Number 
for the Hn region in the catalogue of each galaxy. Column 2: Logarithmic Ha luminosity. Column 3: Shell radius. Column 4: 
Measured shell expansion velocity. Column 5: Shell electron density. Column 6: Ionized mass of the shell. Column 7: Kinetic 
energy of the shell. Column 8: Combined kinetic energy of OB stellar winds, assuming an equivalent number of 03(V) type 
stars. 
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Fig. 1. Left: Emission line profile in Ha of the H n region 403 of M 5 1 in Rand's photometric catalogue. The profile was obtained 
from the Fabry-Perot data cube of the H n region. The Gaussian decomposition shown in the figure was obtained using the task 
PROFIT in the GIPSY package (van der Hulst et al. 1992). Right: Emission line profile of the H n region 3 12 of M 5 1 in Rand's 
catalogue obtained and fitted in the same way as in the left figure. 



stellar winds and also by the supernovae integrated over a time 
scale of 1.0xl0 7 yr, and have studied the evolution of a shell 
over this time scale. Estimates of the ages of H n regions based 
on predicted ratios of emission line strength to continuum level 
(Copetti et al. 1985; Bresolin & Kennicutt 1997) yield ages for 
H ii regions in the range from 2Myr to 8Myr, and show that the 
dynamical timescale adopted here is reasonable. 



3.2. Two basic driving mechanisms: Stellar Winds and 
Supernova explosions 

There are two basic mechanisms which can drive the forma- 
tion of a shell expanding into the ISM from an OB association: 
the stellar winds from young stars and the explosions of super- 
novae. The first mechanism was modelled when accounting for 
the shells observed in the H n region populations of the three 
spiral galaxies described in Paper I. We compared there the 
wind energy input from the equivalent numbers of 03(V) stars, 
derived using the Ha luminosity of a region, with the shell 
kinetic energy inferred from the kinematic observations. The 
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conclusion was reached that shells with radii of up to ~ 0.2 R les , 
(where R lee is the radius of the H n region) could be quantita- 
tively explained by the interaction of the stellar winds with the 
ISM, while some of the bigger shells located within the most 
luminous H n regions show kinetic energies which are too high 
to be explained as the result solely of wind interaction. 

We have here gone on to improve the method used in Paper 
I, and now use the input energy from the stellar winds (here 
derived from the models implemented in Starburst99), but also 
add in the energy released in the subsequent supernovae. We 
are thus modelling the situation in which both mechanisms op- 
erate on the ISM surrounding a star cluster. The winds provide 
a continuous energy source which lasts for the stellar lifetime, 
i.e. ~10 6 yr, by this time an expanding shell has formed within 
the H ii region. This is a single shell even if the stars "switch 
on" at dispersed times, because if the first star causes a shell to 
form, the wind from a second star will expand more quickly un- 
til it reaches the shell, and then couple its impulse to the shell, 
and this will occur for each subsequent star. Then, the super- 
novae occur inside the partially attenuated bubble produced by 
the winds. The situation has been described in some detail by 
Bruhweiler et al. (1980) and by Dyson (1981); the shock pro- 
duced by the SN explosion sweeps up the material inside the 
cavity formed by the stellar winds, and impacts directly on the 
shell; then the strength of the transmitted shock is reduced by 
a factor which depends on the contrast in density between the 
cavity and the newly formed shell. The radiative cooling time 
within the shell can be very short, and then the momentum con- 
servation is applied as the best approximation to describe the 
dynamics of the new formed shell. 

The initial velocity is taken directly from the Ha observa- 
tions, so that Vo is the velocity of the shell as observed within 
an H ii region, in Ha (column 4 of Table 1). The radius at which 
the SNe explode (Ro), is taken to be the radius of a wind blown 
shell. This is estimated applying the model of Dyson (1980) for 
the adiabatic expansion of a stellar wind bubble. The inputs re- 
quired in this model are the energy flux in the stellar winds, the 
initial gas density, and the time during which the stellar winds 
act. Taking the stellar wind energy from Starburst99 models we 
apply the equations of the adiabatic model to a range of values 
of n and to. The optimum combination of these variables for 
producing the measured parameters of the wind blown portion 
of the shell expansion is n=0.5cirT 3 and to = 10 6 yr. Using 
these values the shell radius is found to be ~ 0.3R reg and the 
velocities just slightly lower than the observed Ha velocities. 
In the end we chose to use Ro = 0.5R reg for the initial ra- 
dius of the shell in the post SN phase of our modelling. This 
is the shell radius at which the SNe explode and although it is 
a little higher than the radius derived in the adiabatic model, 
it allows us to simplify our modelling by separating the adia- 
batic expansion phase due to the winds from the momentum 
conserving phase after the supernovae have exploded. We used 
an ambient interstellar density of n=0.1cirT 3 , slightly lower 
than the value estimated from the adiabatic model. This val- 
ues is chosen as a practical compromise which takes partially 
into account that the density decreases with distance from the 
plane of the galaxy, a circumstance of relevance especially 
for the larger HI shells. We confirmed that using a density of 



0.5 cm 3 changes the magnitude of the results by a relative 
modest fraction (slightly (~20%) higher radii and higher ve- 
locities (-13%)). 

Using the model assumptions described above we start with 
the shell whose properties we have derived from our Ha ob- 
servations and run the models adding the supernova inputs as 
explained. We assume that the stellar winds have acted over a 
time interval of t Q = 10 6 yr and that at this time the supernovae 
explode with no time intervals between them. The supernovae 
then sweep up the material within the bubble and push out- 
wards the shell created by the stellar winds. The application of 
momentum conservation gives the radius and the shell veloc- 
ity as a function of time according to the expressions (Dyson 
1980; Bruhweiler et al. 1980): 



Rshell(t) = [Ro + 



^ (t _ to)] (l/4) 

nmeTr 



Vshell(t) = R 



3 3MqVq 

shell 4nm H 7r 



(1) 

(2) 



where Ro and Vo are the radius and velocity of the expanding 
shell at the switch-on time of the supernovae, n is the ambient 
interstellar density and Mo is the mass swept up by the effect 
of the supernova explosions. In Fig. [2] we show the time evolu- 
tion of the shell radius and shell velocity for the H n region 8 
in NGC 1530 (first raw of Table [B- The figure, that can be 
taken as a characteristic example, shows qualitatively how the 
shell forms in the wind blown phase of the expansion. It then 
receives a boost from the supernova input which initially accel- 
erates it, but the velocity is subsequently reduced as the shell 
picks up mass from the surrounding ISM and expands reaching 
a radius within the range measured for HI supershells. 

The mass Mo in Eq. [2] is the total mass contained within 
the bubble at the initial radius Ro, which is assumed to be com- 
pletely swept up as a consequence of the SN explosions. To 
obtain the mass within the bubble we have used the method to 
derive the mass of gas in an H n region described in Relano et 
al. (2005b), where we included a factor to take into account 
the neutral gas not observed directly in Ha (see Giammanco 
et al. 2004). We then apply Eqs. 1 and 2 and compute a "fi- 
nal" radius for the shell, R s h e ii, and a "final" velocity, v s h e ii, at 
a time of 1.0 x 10 7 yr, which agrees with the dynamical age 
range of the suppershells observed in HI (McClure-Griffiths et 
al. 2002; Walter & Brinks 1999). An integration of the models 
at a time of 0.5 x 10 7 yr will produce velocity values between 
15-30kms _1 , while at longer times (5.0 x 10 7 yr) the velocities 
will be ~3kms~ 1 , much lower than those shown here. These 
last values are not expected to have a clearly defined phys- 
ical meaning, since the shell loses its identity at v~5kms~' 
(Dyson 1981), and for the corresponding radii, of order 500 pc 
of greater, a shell will tend to be disrupted by global galactic ef- 
fects such as differential rotation. We present the results of the 
models in Table |2j where we also show the total combined in- 
put energies of stellar winds and supernovae from Starburst99, 
and the kinetic energy of the bubble obtained using the standard 
expression = jm H iV^ r An interesting supplementary exer- 
cise was to run models in which a second burst of star forma- 
tion boosted the injection of energy into the expanding shells. 
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01234 01234 
time (Myr) time (Myr) 

Fig. 2. Left. Evolution of the shell radius of the adiabatic and momentum conserving phases for the parameters shown by H n 
region 8 of NGC 1530. RAd is the initial radius found from the adiabatic model, which predicts in general shell radii of 0.3R reg , and 
velocities comparable with those observed for the ionized shells within the H n regions. Ro is the input radius for the momentum 
conserving phase. We have not explicitly modelled the transition between the phases, and this is represented schematically here 
by the line joining the two curves in the time interval between 1 and 1.2 Myr. Right. Evolution of the shell velocity for the 
same example H n region. The shell velocity predicted in the adiabatic phase is close to the observed velocity of the Ha shell 
(see column 4 of Table [TJ. In our simplified model of the momentum conserving phase, in which the SN input is instantaneous, 
the initial velocity would be boosted to the unrealistically high value indicated in the figure. A schematic representation of the 
transition is given by the line joining the two curves in the time interval between 1 and 1.2 Myr. Although the assumption of 
instantaneous energy injection is too simple, the overall energetics of the model, which predicts a final HI radius of 653.5 pc and 
a corresponding HI expansion velocity of 17 km s (see Tabled are physically reasonable 



According to the time at which the model output is examined 
this produces HI shells either with a range of radii similar to 
those observed but with a velocity range substantially higher, 
or with a range of velocities similar to those observed but with 
a range of radii substantially higher. This suggests strongly that 
the stars in OB associations are coeval. In the next section we 
compare our results with a sample of HI shells observed in two 
different galaxies. 

3.3. Comparison with HI observations 

We have selected a representative sample of HI shells with dif- 
ferent sizes observed by Walter & Brinks (1999) in the dwarf 
galaxy IC 2574 and by McClure-Griffiths et al. (2002) in the 
Milky Way. We have limited our sample to those shells classi- 
fied by Walter & Brinks as "classic holes" and we have avoided 
"merged shells" and chimneys in the McClure-Griffiths et al. 
data set. Our results in Table [2] show good broad agreement 
with the parameters of the observed HI shells reported in 
the two articles cited. The modelled range of shell velocities, 
v s heii = (8 - 17)kms _1 , falls nicely within the range of the 
observed sample, v s h e ii = (4 - 21)kms~ 1 . The range of shell 
masses produced in the models, M s h e ii = (3 - 29) x 10 5 M o is 
within the range of masses observed by McClure-Griffiths et 
al.(2002), (1 - 44) x 10 5 M o but is somewhat higher than the 
range found by Walter & Brinks (1999), (0.02 - 18) x 10 5 M o . 
The first study yields values for the HI number density 0.5- 
1.7 cm 3 , while the second study takes the scale height of 
the galaxy disk into account (the galaxy is in any case much 



smaller) and finds values in the range 0.05-0. 17 cm 3 . The 
smaller values for the range of the basic parameters in Walter 
& Brinks (1999) can be ascribed to the fact that they were 
observing a dwarf galaxy where the range of ISM densities 
is lower, and this will necessarily be reflected in the ranges 
of shell mass and kinetic energy. Our H n region observations 
have all been taken from regions in essentially Milky Way mass 
galaxies. The better agreement between the range of parameters 
derived by applying dynamical modelling to shells observed in 
these regions, using the appropriate values of mass, velocity 
and ISM density, and the HI shells observed in the Milky Way 
by McClure-Griffiths et al. is not fortuitous. 



4. Final Remarks 

It is not surprising to find expanding shells around OB associ- 
ations, but a variety of scenarios have been proposed in which 
these can arise. OB stars have a major dynamical impact on 
their surroundings, which can take four forms: the direct influ- 
ence of ionizing radiation on the surrounding gas, the influence 
of non-ionizing radiation on the surrounding gas, mediated by 
its dust content, the effect of the outflow of gas from the stars 
themselves in the form of stellar winds, and the effects of the 
violent outflow of gas in supernova explosions. Our underly- 
ing concepts here are based on the seminal article by Dyson 
(1981), who considered the dynamical effects of stellar winds 
on H ii regions in the necessary presence of the ionizing radia- 
tion which forms them. 
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Table 2. Energetics of the HI predicted shells for the sample of H n regions in NGC 1530, NGC 3359, NGC 695 1 and NGC 5 194. 
The stellar winds act over 10 6 yr, the time at which the SNe explode, the final integration time of the models is 10 7 yr. Column 1: 
H ii region number. Column 2: Shell radius. Column 3: Shell expansion velocity. Column 4: Mass of the shell. Column 5: Energy 
input from the supernovae and stellar winds. Column 6: Kinetic energy of the shell. Column7: Efficiency of the SN to convert 
their energies into kinetic energy in the ISM. 



Region 


Rshell 


Vshell 




E(SW+SN) 


E K (shell) 


Efficiency 


(number) 


(pc) 


(kms -1 ) 


(10 S M G ) 


(10 52 erg) 


(10 51 erg) 


(%) 


1530-8 


653.5 


17.3 


28.9 


36.1 


8.6 


3.7 


1530-22 


505.8 


13.4 


13.4 


16.8 


2.4 


2.2 


1530-92 


351.6 


9.4 


4.5 


2.7 


0.4 


2.3 


3359-6 


550.2 


14.6 


17.2 


18.1 


3.7 


3.2 


3359-42 


353.5 


9.3 


4.6 


4.7 


0.4 


1.3 


3359-92 


299.7 


8.0 


2.8 


2.0 


0.2 


1.4 


6951-2 


579.4 


15.4 


20.1 


24.1 


4.8 


3.1 


6951-18 


386.8 


10.3 


6.0 


6.0 


0.6 


1.6 


6951-41 


326.9 


8.7 


3.6 


2.7 


0.3 


1.6 


5194-403 


449.7 


11.9 


9.4 


10.5 


1.3 


2.0 


5194-312 


400.9 


10.7 


6.7 


5.5 


0.8 


2.1 


5194-416 


348.8 


9.2 


4.4 


3.7 


0.4 


1.5 



In treating the development of the wind-blown shells we 
must then take into account that on a well determined timescale 
there will be an additional dynamical input from the super- 
novae in the OB association. The timescale is of key impor- 
tance, because during the main sequence lifetime of an O star, 
of order 10 6 yr, a wind-blown shell expanding at a characteris- 
tic velocity of ~50kms~ 1 will reach a radius of ~50 pc. This 
implies that the supernovae, which explode in the star cluster 
will at first expand very rapidly into a low density environment 
cleared by the winds, but will then impinge on the denser region 
marked by the shell and the surrounding ISM. Under these cir- 
cumstances the adiabatic phase (Sedov 1959; Chevalier 1977) 
of the SN remnant expansion will immediately give way to the 
snowplough phase (Chevalier 1974) and the dynamical cou- 
pling of the SN remnant to the expanding shell will be deter- 
mined by momentum conservation. We find that a typical frac- 
tion of the expansion energy of the SN remnant transferred to 
the shell will be (1-4)% (see last column in Tabled, with the 
rest either taken up as turbulence or radiated away. Further ob- 
servations and more complete analysis would be necessary in 
order to set more rigorous constraints on the efficiencies under 
a variety of conditions. For example, if the largest bubbles blow 
out when expanding into lower gas densities, the efficiency of 
the process will be predictably lower than those found here. 

To summarize, the detection and measurement of the phys- 
ical parameters of dominant expanding shells within a major 
fraction of highly luminous extragalactic giant H n regions pre- 
sented in Paper I and fully confirmed by the observational re- 
sults in the present article, has given us a starting point from 
which to present a scenario relating the presence of these shells 
to the presence and the properties of the much larger HI shells 
detected in the Milky Way and in nearby galaxies. The result 
of applying a straightforward model embodying this scenario 
is that the HI shells in general most probably had Ha shells as 
their precursors. The time scales within the model are such that 
one should not expect to observe both the precursor shell and 
the HI shell for the same system, so that our results are neces- 



sarily statistical in nature. However we can conclude that the 
observations of Hatzidimitrou et al. (2005), who observed stars 
in OB associations at the centres of only a minority of their 
HI shells, are fully compatible with our scenario. A suggestive, 
though perhaps not conclusive secondary inference, is that the 
stars in an OB association probably form in a single group not 
dispersed in formation epoch. 
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